Modified holographic Ricci dark energy coupled to interacting dark matter 
and a non interacting baryonic component 
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We examine a Friedmann- Robertson- Walker universe filled with interacting dark matter, modified 
holographic Ricci dark energy (MHRDE), and a decoupled baryonic component. The estimations of 
the cosmic parameters with Hubble data lead to an age of the universe of 13.17 Gyr and show 
that the MHRDE is free from the cosmic-age problem at low redshift (0 < z < 2) in contrast 
to holographic Ricci dark energy (HRDE) case. We constrain the parameters with the Union2 
data set and contrast with the Hubble data. We also study the behavior of dark energy at early 
times by taking into account the severe bounds found at recombination era and/or at big bang 
nucleosynthesis. The inclusion of a non interacting baryonic matter forces that the amount of dark 
energy at z t ~ 0(1) changes abruptly implying that Q x (z ~ 1100) = 0.03, so the bounds reported 
by the forecast of Planck and CMBPol experiments are more favored for the MHRDE model than 
in the case of HRDE cutoff. For the former model, we also obtain that at high redshift the fraction 
of dark energy varies from 0.006 to 0.002, then the amount of Q x at the big bang nucleosynthesis 
era does not disturb the observed Helium abundance in the universe provided that the bound 
Q x (z ~ 10 10 ) < 0.21 is hold. 
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I. INTRODUCTION 

The holographic principle states that the maximum 
number of degrees of freedom in a volume should be 
proportional to the surface area pj, 0], [||, Q. Using 
the effective quantum field theory it was shown that the 
zero-point energy of a system with size L should no ex- 
ceed the mass of a black hole with the same size, thus 
L 3 p\ < LMp, where p\ corresponds to the quantum 
zero-point energy density [f| and Mp 2 = 8nG. The lat- 
ter relation establishes a link between the ultraviolet cut- 
off, define through p\, and the infrared cutoff which is 
encoded by the scale L. Applying this novel principle 
within the cosmological context implies that the dark en- 
ergy density of the universe p x takes the same form of the 
vacuum energy, p\ — p x . Using the largest L as the one 
saturating the above inequality, it turns out to be the 
holographic dark energy is given by p x = 3c 2 M 2 p L~ 2 , 
where c is a numerical factor. The IR cutoff has been 
taken as the large scale of the universe, Hubble horizon 
@ , particle horizon, event horizon 01 or generalized 
IR cutoff @, ||, 0, 0, 0, 01, 0] amongst 
many others. One of the main reasons in working within 
the framework of dynamical dark energy such as a HDE 
relies on the need of explaining the current accelerated 



a Electronic address: chimcnto@df.uba.ar 
b Electronic address: monicafortc@fibcrtel.com.ar 
c Elcctronic address: martin@df.uba.ar 



phase of the Universe; this fact has been confirmed by a 
plethora of observational tests such as high redshift Hub- 
ble diagram of type la supcrnovae as standard candles 
00 and accurate measurements of cosmic microwave 
background (CMB) anisotropies pol - E^ . In order to play 
the role of a dynamical dark energy model, the infrared 
cutoff will be considered as a function of the cosmic time 
so the holographic dark energy will evolve dynamically. 
Here, we will focus our attention on an extended version 
of the well known Ricci scalar cutoff 0] . An important 
feature of this model refers to the fine-tuning problem, 
that is, due to the dark energy density is based on space 
time scalar curvature, without involving a Planck or high 
physical energy scale, the fine-tuning problem is avoided 
and the coincidence problem is also discarded within this 
context Several works have been devoted to obtain 
cosmological constraints on the holographic Ricci dark 
energy model 0] , [24| , Phi . , 0] or generalized ver- 
sions of the latter one [2q |. Nevertheless dark energy 
is not the only mysterious element in the Universe, the 
necessity of a dark matter component comes from as- 
trophysical evidences of colliding galaxies, gravitational 
lensing of mass distribution or power spectrum of clus- 
tered matter 0], 0. Moreover, the astrophysical ob- 
servations from the galactic to the cosmological scales 
indicate that dark matter is a substantial component to 
the universe's total matter density, being responsible for 
the structure formation in the Universe [27j . Presents at- 
tempts to understand the physics behind the dark sector 
composed of dark matter and dark energy have indicated 
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that there is an unavoidable degeneracy between dark 
matter and dark energy within Einstein's gravity, namely, 
there could be a hidden non-gravitational coupling be- 
tween them without violating current observational con- 
straints and thus it would be interesting to develop ways 
of testing exchange of energy in the dark sector. More 
precisely, if dark energy interacts with dark matter, there 
is a change in the background evolution of the Universe 
that allows us to constrain a phenomenological type of 
interaction. Therefore, an holographic scenario becomes 
a physically viable model when one takes into account a 
possible interaction between the dark matter and dark 
energy. It entails that the dark matter feels the presence 
of the dark energy through the gravitational expansion of 
the universe plus the exchange of energy between them. 
In fact, we will follow a phenomenological approach by 
studying the properties hidden in a particular kind of in- 
teraction and then one confronts the theoretical model 
with the available observational data. Recently several 
known linear and nonlinear interactions in the dark sector 
have been generalized [2^| . It was introduced an effective 
one-fluid description of the dark components and shown 
that interacting and unified models arc related to each 
other. It should be stressed that interacting dark e nerg y 
scenarios have been studied by many authors [2!| , [30j . 

It is well known that some new physics may be show- 
ing up at high rcdshift taking into account strict limits 
coming from big bang nucleosynthesis (BBN) data. More 
precisely, the 4 i7e abundance has often been used as a 
sensitive probe of new physics. This is essentially due to 
the fact that nearly all available neutrons at the time of 
BBN, in a scale of lMev of temperature, at the time of 
BBN end up in 4 He and the neutron-to-proton ratio is 
very sensitive to the competition between the weak inter- 
action rate and the expansion rate. For example, a bound 
on the number of rclativistic degrees of freedom (d.o.f), 
<?* , at the time of BBN commonly known as the limit on 
neutrino flavors, N v , is derived through its effect on the 
expansion rate, H ex ^fglT 2 where <?* = 2 + 7/2 + 7N u /4: 
counts the relativistic dof in photons, e ± pairs, and N v 
neutrino species; having assumed a Universe dominated 
by radiation so that p cx <7*T 4 [3l|. Moreover, the pres- 
ence of vacuum energy during BBN is well motivated 
both by considerations of dark energy as well as infla- 
tion, giving as stringent bound O x (lMcv) < 0.21 [3lj ] - 
Besides, the physics at recombination era gives also some 
constraints in the amount of dark energy at such primor- 
dial era which also has to be consistent with the severe 
bounds provided by BBN data mentioned above. In par- 
ticular, stringent signal could arise from the early dark 
energy (EDE) models, that is, uncovering the nature of 
dark energy as well as their properties to high redshift 
along with their effects imprinted on the universe could 
provide invaluable guide to the physics behind the re- 
cent speed up of the universe HH • Therefore, any serious 
dark energy model used for constraining the present-day 
value of fl x (z = 0) also has to be consistent with the 
early bounds on the fraction of dark energy at primordial 



eras such as recombination or BBN. Not too long ago it 
was examined the current and future data for constrain- 
ing the amount of EDE, the cosmological data analyzed 
has led to an upper bound of Cl x (z ~ 1100) < 0.043 
with 95% confidence level (C.L.) in case of relativistic 
EDE while for a quintessence type of EDE has given 
Q x (z ~ 1100) < 0.024 although the EDE component 
is not preferred, it is also not excluded from the current 
data [32|. Another appealing forecast for the bounds of 
the EDE taking into account Planck and CMBPol ex- 
periments can be found in (33|. More precisely, assum- 
ing a fl x (a ~ 10 -3 ) ~ 0.03 among other priors, it was 
checked the stability of these values, interesting enough 
was the la error coming from Planck experiment giving 
as result a ^ lanck ~ 0.004 whereas the CMBPol improved 
this bound by a factor 4 „CMBPoi _ 1Q -3 Besides, 
some new limits on EDE from the CMB using the data 
from the WMAP satellite on large angular scale and the 
South Pole Telescope (SPT) on small angular scale were 
obtained in [34|. Considering the CMB data alone it got 
a bound of Cl x (z ~ 1100) < 0.0018 very similar to the one 
reported in [33j . In addition, the constraints on the vari- 
ation in the fine structure constant [35[ in the presence 
of EDE gave an upper bound of Vl x (z ~ 1100) < 0.06 at 
95% C.L. which is weaker than the bounds reported in 

[H-m. 

The present article is outlined as follows. We inves- 
tigate a universe composed of interacting dark matter, 
modified holographic Ricci dark energy (MHRDE), where 
the exchange of energy in the dark side is proportional to 
derivative of total dark sector energy density, and a de- 
coupled component that could behave as baryonic matter 
at early times. We use the Hubble data and the Union2 
compilation of SNe la for constraining the cosmological 
parameters, thus we compare phenomenological aspects 
of MHRDE and HRDE holographic dark energy models. 
We also make a kinematic analysis for studying the be- 
havior of decelerating parameter, equations of state and 
the ratio dark matter to dark energy. Using the best fit 
values, we estimate the age of the universe and explore 
the cosmic age-problem. As a complementary tool for 
getting more accurate constraints on both models, we 
examine the behavior of dark energy at early times. 

II. THE MODEL 

We assume a flat FRW universe filled with three dif- 
ferent components, an interacting dark sector composed 
of a nearly pressurelcss dark matter, MHRDE, and a de- 
coupled baryonic contribution with energy densities p Cl 
p x , and pb, respectively. We adopt as equations of state 
w c = Pc/pc for dark matter, u> x = p x /p x for dark energy, 
and uib = Pb I Pb for radiation, thus the Einstein equations 
read 

3H 2 = p c + p x + p b , (1) 
p c + p x + 3H[(u> c + l)p c + (uj x + l)p x \ = 0, (2) 
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pb + 3H(uj b + l)p b = 0, 



(3) 



where a is the scale factor, H = a/a stands for the Hub- 
ble expansion rate. Here, we will use the holographic 
principle within the cosmological context by associating 
the infrared cutoff L with the dark energy density, thus 
we take L~ 2 in the form of a linear combination of H and 
H 2 0: 



Px = 



a- P 



H 



3a 



H- 



(4) 



being a and P two free constants. In particular, we obtain 
p x oc R for a = 4/3 [lj], where R = 6(H + 2H 2 ) is the 
Ricci scalar curvature for a spatially flat FRW space- 
time. 

The use of the variable t] = ln(a/ao) 3 , where ao is set 
as the value of the scale factor at present, allows us to 
rewrite Eqs. ©-Q as 



P = Pc 



Px 



-{u> c + l)p c - (U) x + l)p x , 



p = -ap c 



Ppx, 



Pb = -(u>b + l)Pb, 



(5) 



(6) 



(7) 



(8) 



where the prime stands for derivatives with respect to 
the new variable ' = d/drj and the condition < j3 < a 
is imposed to avoid a phantom scenario. From Eq. ((SJ) 
is clear that the radiation component is decoupled from 
interacting dark sector, so the exchange of energy only 
takes place between the dark matter and dark energy, 
thus pb = pboa~ 3 ( u ' b+1 ' > and its density parameter is f2& = 
Pb a- 3 ^ +i y(3H 2 ). 

MHRDE (U), with a term proportional to H leads 
to Eq. 0, which looks like a "conservation equation" 
for the two dark components with constant coefficients. 
We will refer to the Eq. ([7]) as the modified conser- 
vation equation (MCE). Comparing the whole conser- 
vation equation (WCE) flHJ} and the MCE 0, namely 
(u> c + l)p c + (w s + 1)/5 X = a Pc + fipxi we obtain the com- 
patibility relation 



(a-uj c -l)r + /3-l, 



(9) 



between the equation of state of both components and 
its ratio r = p c /p x . In what follows, we will use the 
MCE ([7]) with constant coefficients a and (3 instead of 
the WCE © with non-constant coefficients. In some 
sense, the WCE @ and the MCE give rise to different 
representations of the mixture of two interacting dark 
fluids and clearly these descriptions are related between 
them by the compatibility relation ©. Therefore, the 
MHRDE conveniently links a model of two interacting 
fluids having variable equations of state with a model 



of two interacting fluids with "constant equations of 
state" . 

Using Eqs. (J5J5]) the total pressure is p = p c + p x and 
the effective equation of state of the dark sector (EOS), 
uj = p/p can be rewritten as 



-P-P 



uj c r + uj x 
1 + r 



(10) 



At this point, we introduce an interaction ZHQ be- 
tween the dark components by splitting the MCE ([7]) in 
the following way 



ape 



-Q, p' x + p Px = Q. 



(11) 



Now, we assume a pressureless dark matter [ui c = 0), 
hence the equation of state of dark energy © becomes 
linear in r 



(a - l)r + P - 1. 



(12) 



The next step is to introduce a phenomenological inter- 
action between the dark components in order to extract 
some physics information about the behavior of them. 
We are going to study an interacting scenario where the 
exchange of energy between dark matter and dark en- 
ergy is proportional to p' . We will employ the method 
developed by one of the authors [2£| based on the source 
equation for obtaining the total energy density of the 
dark sector once the interaction is given, then we will be 
able to reconstruct the partial energy densities. The new 
kind of interaction was introduced in (29| and reads as 

(u s + l-a)(u s + l-p) , 

Q = / — — tta p> ( 13 ) 

(Us + 1)A 

where u> s is a constant that varies between a — 1 and 
p — 1, A = a — P, and Q < 0. Taking into account 
that the partial energy densities p c and p x appears as a 
linear combination in the conservation equationQ, the 
interaction (|13|) can be expressed as a linear combination 
of p' c and p' x also. In what follows we will employ the 
method of the "source equation" developed in [29| for 
obtaining the total energy density as 



bia Vs+i) 



. ko-30-.+i) 



(14) 



where b\, 62 are integration constants. At early times 
the effective energy density of the dark sector takes the 
form p ~ bia~ 3al3 /( Us+1 \ in order to have a term a -3 in 
the dark sector we need to take a/3 = u> s + 1. In latter 
case, the dark matter and dark density parameters take 
the forms 



n r . = 



(l - P)h + (u s - P + l)b 2 a 
{b 1 + b 2 a- 3 ">)A 



— 3u) s 



l)bi 



-3oj s 



(h +b 2 a~ 3 "°)A 



(15) 



(16) 
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Dark matter and dark energy densities behave as a~ 3 
with a constant ratio r e ~ (1 — f3)/(a— 1) at early times. 
However, at late times the dark components behave as 



-3a/3 



so the parameter densities ([15]) and (fTB|) give Q c 



[uj s -fi + I] /[a - 13] and Q x ~ [a - 
n ~ /3(a - l)/a(l - (3) and w x 
l) 2 ]/a(l-/3). 



1 — cj s ]/[a — hence 
~ [/3(a - l) 2 - a((3 - 



III. COSMOLOGICAL CONSTRAINTS 

In what follows, we will place some constraints on the 
model, mentioned in the last section, using the observa- 
tional Hubble H(z) data and the constraints imposed by 
the Union 2 compilation of SNe la. The function H{z) 
plays a crucial role to understand the properties of the 
dark energy since its value is directly obtained from astro- 
physical observations. More precisely, the differential age 
data of astrophysical objects that have evolved passiv- 
elly during the history of the universe (e.g. red galaxies) 
allows to test theoretical cosmological models through 
the predicted Hubble function H(z) = —(1 + z)~ 1 dz/dt 
expressed in terms of the redshift z. Hence, we obtain 
the function H(z) by direct determination of dz/dt [36j . 
This can be achieved by identifying some "clock" galax- 
ies that exhibit a uniform distribution of star population 
The 12 observational H(z) data is listed in [37| . 
There, H b s (zi) and H j, s (zk) are uncorrelated because 
they are obtained from the observations of galaxies at dif- 
ferent redshifts, where z varies over the interval [0, 1.75]. 
The statistical analysis is based on the x 2_ function of the 
Hubble data which is constructed as (e.g.[39[) 



[H(9,z k 

fc=i 



H bs(zk)Y 



a(z k Y 



(17) 



where the 9 symbol refers to the set of cosmological pa- 
rameters, H f )s (zi c ) is the observational H(z) data at 
the redshift Zk, a(zk) is the corresponding la uncer- 
tainty, and the summation is over the 12 observational 
H(z) data listed in [37}, [H|. From this quantity, the 
probability distribution function (PDF) is constructed as 



V = We'*'/ 2 where W is a normalization factor. 

The x 2 -function will be minimized for obtaining the 
best-fit values of the random variables 8 C that correspond 
to a maximum of V . The best fit parameters 9 C are those 
values where Xmin(@c) leads to a local minimum of the 
X 2 (^-distribution. If X 2 d . . f = X^MKN -n) < 1 the 
fit is good and the data are consistent with the considered 
model H(z; 0)[39|, where N indicates the number of ob- 
servational data whereas n counts the number of param- 
eters so x 2- fu n ction has N — n degrees of freedom. To 
better understand the cosmological constraints coming 
from x 2_ statistical method employed here, we are going 
to place constraints over all parcmctcrs, taken in pairs, 
while the others are taken as priors, namely, we con- 
sider one pair of unknown parameters and obtain their 
mean value, then we choose another pair of parameters 
and repeat the process until all the set of (9-paramctcrs 
have been properly estimated. From the latter analysis, 
we are going to obtain the 68.3% and 95.4% confidence 
level (C.L.) of probabibility that in the case of two inde- 
pendent parameters corresponds to the random data sets 
which satisfy the inequality Ax 2 = X 2 {@) — Xmini^c) < 
2.30 and Ax 2 < 6.17 respectively, these contours are usu- 
ally closed ellipses. Here, iV = 12, n = 2, the string of 
parameters is 9 = (Hq, a, j3, f2 c o, Q x o, ^t>o), an d the the- 
oretical Hubble function in terms of redshift z is given 
by 

H{z)^H [A(l + z) 3 + B(l + z) 3 ^+ 1 '>+Ci b0 {l + z) :ia } 1 / 2 

(18) 

where the constants A and B are written in term of pa- 
rameters as 



_ {Us + 1 - oQ^cO + (Us + 1 - P)QxO ^ ig ^ 



B 



(a-l)ncO + G8-l)fi*o 



(20) 



having used the standard definition of the density param- 
eters Qio = pio/SHy with i = {x, c, b} and the flatness 
condition f2 c o + Q x q + fibo = 1 is hold. 



Applying x 2 -method to the theoretical Hubble (|T8|) gives us the confidence level associated with the la and 2a 
probability for all the possible pairs (see Fig. [TJ, the best-fit value for each pair is represented by a dot and the 
dashed zones exclude places (in the parameter space ) where the conditions f^o > fail to be guaranteed, the range of 
the parameters a and (3 lead to a phantom scenario or the parameter densities ttio take values that are not consistent 
with the literature (see FigfT]). The best- fit values for each pair of parameters and the corresponding values of 
Xdof P cr degree of freedom are gathered together in Table H] Up to this point, talking in broad terms, we found that 
the values of 9 show a small variation around their mean values: Hq = 73.89lg'g4km s _1 Mpc _1 , a = 1.055io!o38) 
(3 = 0. 126±g:g||, n c0 = 0-18±n;8L tt x0 = 0.768ig;^g and O fc0 = 0.038^;^. The value of H obtained is close to the 
one reported by Riess et al [4y] , Hq = 74.2±3.6kms~ 1 Mpc~ 1 at 68% C.L., being the measurement of Hq obtained from 
the magnitude- redshift relation of 240 low- z Type la supernovae at z < 0.1; the absolute magnitudes of supernovae 
are calibrated using new observations from Hubble Space Telescope (HST) of 240 Cepheid variables in six local Type 
la supernovae host galaxies and the maser galaxy NGC 4258. Further, a 7-year WMAP analysis prefers, but does not 
directly measure, H = 71.0 ± 2.5kms" 1 Mpc- 1 j43|- The value of the Hubble constant was evaluated over a range of 
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FIG. 1: The 68.3% and 95.4% confidence level contours for all pairs of 9 — (Hq, a, /3, Q c , fl x , O^) is depicted. In each figure, the dot indicates the 
best fit observational value obtained with the H(z) function, and the dashed bands denote the forbidden regions for the parameters. 
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(1.017,0.11) 


(1.07.0.18) 
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(0.17.0.77) 


(0.16,0.04) 


(0.17.0.77) 


(0.19.0.03) 


(0.77,0.04) 




0.779 


0.709 


0.707 


0.707 


0.707 


0.705 


0.703 


0.703 


0.702 


0.702 


0.702 


0.702 


0.702 


0.702 


0.702 



TABLE I: The best-fit values for each pair of parameters and the corresponding values of Xd.o.f P er degree of freedom are 
indicated. The average values obtained from these partial adjustments are Ho = 73.89, a = 1.05, f3 = 0.12, £l c o = 0.18, 
0*0 = 0.768 and Q b0 = 0.038 
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redshifts 0.03 < z < 0.5, assuming a cosmological model with f2 c o = 0.27 and Q. x q = 0.73. The value determined for 
the Hubble constant is H = 75.9 ± 3.8kms _1 Mpc _1 [13]. Besides, the value of il c0 is consistent with il c oh 2 = 0.11 
and h = 0.7 (43[. We will use these mean values to analyze the main traits of the model with the third component, in 
particular the issue of early dark energy will be addressed with some detail as a way of further constrain our model 
with the physics behind the primordial eras such as recombi nation or big bang nucleosynthesis. 



2D Confidence level 



Priors 



Best fits 



(n m0 , O^o, 0) = (0.21, 0.75, 0.005) (a, h) = 

(Q m0 , n x0 , a) = (0.21, 0.75, 1.0887) (ft h) = 

(Oco, a, 0) = (0.75, 1.0887, 0.005) (fi m0 , h) 

{Q m0 ,a,f3) = (0.21,1.0887,0.005) (Q x0 ,h) 

(Umo,a,h) = (0.21,1.0887,0.835) (fi*o,,8) 

(fixo, a, h) = (0.75, 1.0887, 0.835) (fi m0 , ft 

(Qmo, ft h) = (0.21, 0.005, 0.835) (Q. x0 , a) 

[a, (3, h) = (1.0887, 0.005, 0.835) (0*0, 0„ 

(O m0 , O l0 , h) = (0.21, 0.75, 0.835) (a, /3) = 



(1.089±g;°|5,0.835±S;SS) 
(0.005tg;S,C 



,0.835±g;SS) 

,0JJ35lg:gg) 
0.835i°;^) 

0.0051°:™) 
,0.005t°;° 7 7?) 
l-089+g:gi) 

04 n 91 +0.76-, 
045 u - zl -0.86^ 

(1.089i;;^,0.005t°;tg) 



= (0.75tH 3 2 , 
= (0.75l°;g«, 

= (0.75185, 
o) = (0.75+° 



TARLE II: Observational bounds for the 2D C.L. obtained in Fig. |2j 
by varying two cosmological parameters. It is reported the best fit val- 
ues of the cosmological parameters with theirs corresponding marginal 
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As is well known, distance indicators can be used for 
confronting distance measurements to the corresponding 
model predictions. One of the most useful ones are those 
objects of known intrinsic luminosity such as standard 
candles, so that the corresponding comoving distance can 
be determined. That way, it is possible to reconstruct the 
Hubble expansion rate by searching this sort of object at 
different redshifts. The most important class of such in- 
dicators is type la supernovac. Then, we would like to 
compare the Hubble data with the Union2 compilation 
of 557 SNe la [4l|. In order to do that, we note that the 
apparent magnitude of a supernova placed at a given red- 
shift z is related to the expansion history of the Universe 
through the distance modulus 

/ , = m-M = 51og^+/x , (21) 
h 

where m and M are the apparent and absolute magni- 
tudes, respectively, ^ = 42.38, h = iJ /100km _1 s _1 , 
and <Il(z) = Hq(1 + z)r(z), being r{z) the comoving dis- 
tance, given for a FRW metric by 



r(z) = 



dz' 



(22) 



To confront the model with supcrnovae data set we con- 
struct the corresponding \ 2 estimator 



fe=i 



cr(-Zfc) 5 



(23) 



where TV = 557 and the cosmological parameters are 
9 = (iJ , a, /3, ri c , fia, fih). Using the Union2 data set, 
we will obtain nine two-dimensional confidence contours 



associated to la and 2a error (see Fig. ^). Thus, 
we obtain the best fit values for nine cases and calcu- 
late the corresponding marginal la error bars [42j as 
it can be seen in Table (JTTJ) . The dashed zones are ex- 
cluded from the analysis due to different reasons such 
as it can be that the range of the parameters a and f3 
lead to a phantom scenario, the parameter densities ftio 
become negative or take values that are not consistent 
with the literature. More precisely, we found that the 
values of a vary over the interval [1.08,1.64] within la 
zone whereas [3 G [0.005; 0.445]. At la C.L. the dark 
energy density parameter at z = goes from 0.34 to 
0.816, fl m0 G [0.21,0.975], and h G [0.832; 0.8350] (see 
Table (|TT| ). The difference between the forecast made with 
the Hubble data and Union2 set is most sharpest in the 
case of f3 parameter, it exhibits a disagree of 0.9%. For h 
and f2 m o the discrepancy between both set does not reach 
0.14% whereas the values of a and fl x o obtained with 
Hubble data disagree with the ones of Union2 by 0.03%. 
In order to corroborate our previous analysis, we also per- 
formed a global statistical analysis with the Union2 data 
set by taking into account a global minimization of the 
five parameters. The latter procedure leads to the best 
fit values (h,a,P,Sl c ,to x ) = (0.83,1.08,0.005,0.21,0.75) 
along with Xd.o.f = < 1j indicating that our esti- 
mations of the cosmological parameters are trustworthy. 
The statistical estimations performed with the Union2 
data set are also consistent with the ones obtained from 
the Hubble data set. 

The density parameters f2 c , £l x , fib, and the ratio r = 
Q c /Q x in term of the redshift z are given by 



_ 1_ / (l-/3)^ + /3(a-l)g(l + z) 3 ^- 1 ) 
c ~ A \A + B(l + zfW-V + n w (l + zf^- 1 ) 

n b0 (i + zf^ 



A + B{1 + z) 3 ^- 1 ) + fi w (l + z) 3(Q - 1} ' 



(24) 
(25) 



1_ / (a- l)A + a(l - /?)ff(l + z) 3 ^- 1 ' 

A \A + B(l + zf^- 1 ) + n w (l + z) 3 («-!) J ' 

(26) 



= (l-p)A + (3{a-l)B(l + zyW-V 

r A + Bii + zfi^- 1 ) + n b0 {i + z) 3{a - 1) 



The aforesaid model (| exhibits dark matter and 
dark energy components with energy densities of simi- 
lar order of magnitude at redshifts z < 2 including the 
present-day scenario (see FigOJ ■ Now, we are interested 
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FIG. 2: Two-dimensional C.L. associated with lcr and 2cr for different 
planes obtained to the Union 2 compilation of SNe la. 



in studying the behavior of kinematic quantities such as 
deceleration parameter and the equations of state. In- 
deed, we readily get to obtain us(z), ui c = 0, ui x (z), 
uif, = a — 1, and q{z) = [1 + 3w(z)]/2 in term of the 
rcdshift 



uj = -1- 



A + apB(l + z) 3 ^- 1 ) + afi M (l + z) 3 ^" 1 ) 



.4 + S(l + z) 3 ^/ 3 - 1 ) + n M (i + z) 3 ^- 1 ) 

(a-l)f! h0 (l + z) 3 



(ff/^o) 2 



(28) 
(29) 



FIG. 3: The density parameters Q c , and SI5 as a function of the 
rcdshift are depicted. Also the ratio Vt c /Vt x is shown. 



It seems like the model experiences a dust-like behavior 
around z = 7.3 and the transition towards the accelerated 
regime takes place at z — 0.87 so it is retarded when the 
nonintcrcating baryonic matter is added. Besides, the ac- 
tual values of all equations of states and the deceleration 
parameter are ojq = —0.63, lu c0 = 0, ui x q = —0.82, uj^ = 
0.06, and qo = —0.45, respectively. These values are con- 
sistent with the high-z supernova data which provide the 
most stringent limit on uj xQ . Using WMAP+BAO+SN, 
it was obtained u x0 = -0.980 ± 0.053 (68% C.L.). The 
error does not include systematic errors in supernovae, 
which are comparable to the statistical error, thus, the 
error in u} x0 from WMAP+BAO+SN is about a half of 
that from WMAP+BAO+ff [ll- The model does not 
cross phantom divide line at any stage of its evolution, 
moreover in the remote future decreases monotonically 
reaching a value of uj x q ~ —0.82, in the same way as 
happens for a quintessence dark energy model, so it does 
not exhibit a quintom behavior (45j . 

The cosmic age-redshift relation for our model reads 



H t(z) 



dx 



{l + x){H{z)/H o y 



(30) 



where the time origin is set at z = oo and the time is 
measured in units of Hq. For the best fit parameters, we 
found that the age of universe is 13.31 Gyr without the 
third fluid whereas its inclusion leads to 13.17 Gyr. Both 
values are very close to the one reported by WMAP-7year 
project, thus it found a 13.75 ± 0.13 Gyr with WMAP 
only and 13.75 ± 0.11 Gyr with WMAP + BAO + H 
[4q . Because the cosmological constraints with the Hub- 
ble data only cover redshifts over the range < z < 2, the 
comparison with cosmic milestones will be trustworthy 
in this range only, for the latter reason we consider two 
old stellar sources such as the 4 Gyr old galaxy LBDS 
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FIG. 4: It shows uj(z), uj c (z), u)j;(z), Ub(z) 7 and q{z) in term of the 
rcdshift z. 



53W069 at redshift z = 1.43 [47[ and the 3.5 Gyr old 
galaxy LBDS 53W091 at redshift z = 1.55 [see Fig[5]. 
We depict the age-redshift relations at the best-fit value 
corresponding to the two models mentioned before (see 
Fig J5|) . We find that the Ricci-like holographic dark en- 
ergy model cannot be accomodated well under the age- 
redshift curves exhibiting a cosmic-age problem at low 
rcdshift, namely, the universe cannot be younger than 
its constituents. On the other hand, the MHRDE seems 
to be free from the cosmic-age problem at low redshift. 

Now, we would like to attend an appealing discussion 
concerning the behavior of dark energy at early times 
within the framework of interacting dark sector for the 
two cases; one case corresponds to the interacting dark 
sector only, whereas the second one is the model with 
the three components. In doing so, we exhibit the cos- 
mological evolution of the density parameters fl x (z) for 
the first model when the universe is filled with dark mat- 
ter and dark energy only, whilst the second model refers 
to a Universe filled with an interacting dark sector plus 
a non interacting baryonic fluid [see Eq. (|26|)]. The dark 
energy is depicted over the range z £ [0, 10 15 ] in order to 
have in mind both the bounds coming from a recombi- 
nation era as well as those produced by BBN data. The 
idea is to compare the behavior of dark energy at early 
times for the two models mentioned in this article, and 
thus, we also will contrast our finding with the ones re- 
ported in the literature. Although both models have the 
same kind of interaction it turned out that the behav- 
ior of their density parameters are very different at early 
times (see Fig|5]). In broad terms, both models exhibit a 
stable behavior about the current value 0.77 at low red- 
shifts ( z < 0.1) and then they begin to separate abruptly 



FIG. 5: It shows the cosmic age-redshift relation for the 
MHRDE with and without the noninteraction baryonic mat- 
ter ( dashed green lines), HRDE (dashed black curve). 



within the interval 0.1 < z < 15, showing a descending 
slope. If the non-interacting fluid is excluded, the den- 
sity parameter Q x exhibits another plateau very similar 
to one of the late times but now the amount of dark en- 
ergy is fixed around 0.23, so it would not strictly satisfy 
the recombination bound, nevertheless at z ~ 10 8 cor- 
responding to a BBN era it would be close to meet the 
condition Q x (z ~ 10 10 ) < 0.21. On the other hand, the 
inclusion of baryonic non interacting fluid makes possi- 
ble that the amount of dark energy continues descending 
rapidly, exhibiting values which are perfectly in agree- 
ment with the stringent constrains provided by the EDE 
at recombination era. For instance Cl x (z = 10 3 ) = 0.038 
and n x (z ~ 1100) = 0.032, these values indicate that the 
second model is consistent with the forecast of Planck 
and CMBPol experiments (33|, meeting below the up- 
per bound provided by the constraints on the variation 
in the fine structure constant [35| and showing a slightly 
discrepancy of one order of magintude with the bound re- 
ported in [34| when the CMB data alone is used. As we 
have already mentioned, the presence of dark energy at 
BBN era should not disturb the observed Helium abun- 
dance in the universe which is regarded as one of the 
major evidence in supporting the big bang theory. In re- 
lation with that, we found that the amount of dark energy 
is il x (z ~ 10 10 ) < 0.21 at BBN [H|. In our model, we ob- 
tained fl x (z = 10 9 ) = 0.0068, and Q x (z = 10 12 ) = 0.0021 
so the fraction of dark energy at early times clearly ful- 
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fills the aforesaid constraint (see Figj6]). It should be 
also stressed that the most importants changes are pro- 
duced near the transition era when the universe entered 
in the accelerated regime (Fig(5J). In addition, we find 
that within the framework of holographic Ricci dark en- 
ergy model the bound at recombination era is not sat- 
isfied with the baryonic component or without it. In a 
future research, we will explore this kind of interaction 
by taking a radiation or baryonic term coupled to the 
dark sector; we will examine the changes introduced in 
the behavior of dark energy at early times [49| . 

Finally, notice that the value of the cosmological pa- 
rameters used here are not arbitrary because these pa- 
rameters turned to be consistent with three important 
data set: i— the present-day scenario constrained with 
the Hubble data and SNe la data set, ii— the recombina- 
tion bounds for EDE, and iii— the BBN data. Another 
useful observational constraints can be found in the last 
scattering surface (LSS), namely, during the galaxy for- 
mation era (1 < z < 3) dark energy density have to be 
subdominant to matter density so it should satisfy that 
< Sl x < 0.5 as happens in our case. 



IV. CONCLUSIONS 

We have examined a FRW universe filled with dark 
matter, a MHRDE with a cutoff given by L~ 2 = (H + 
^■H 2 )/(a — f3) that includes Ricci ansatz for a = 4/3, 
and a baryonic non interacting component decoupled 
from the dynamics of dark sector. We have studied the 
case in which the interaction in the dark sector is pro- 
portional to its energy density derivative, namely, Q oc p' 
and found that the total pressure of both dark compo- 
nents becomes strongly negative, violating the strong en- 
ergy condition in the present epoch. 

We have performed a x 2 -analysis using the Hub- 
ble data and built the corresponding la and 2a C.L. 
(see FigdJ) for each pair of parameters (see Table fl}. 
The mean values of 9 C — (Hq, a, /3, £l x o, S^cOj ^6o) are 



Ho = 73.89i£<gkms 




FIG. 6: The cosmic behavior of the dark energy density £l x in 
term of the redshift over the interval [0, 10 15 ] for the MHRDE 
model. The green solid line shows the case where the universe 
is filled with dark matter and dark energy, whereas the red 
solid line indicates a universe filled with three components, 
two of them encoded in the interacting dark sector, and the 
third one is a baryonic non interacting fluid. 
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along with & Xd o f < 1 • Taking into 



account these 9 C into the density parameters of the dark 
components, we have found that they have similar be- 
havior for redshifts z < 2 including the present-day sce- 
nario (see Figj3j) . Regarding the statistical estimations 
made with the Union 2 compilation of SNe la, we have 
found that a varies over the interval [1.08,1.64] within 
la region whereas /3 G [0.005; 0.445]. At la C.L. the 
dark energy density parameter at z = goes from 0.34 
to 0.816, Q m0 G [0.21,0.975], and h G [0.832; 0.8350] 
(see Fig. © and Table (|H| )). We have also performed 
a global statistical analysis with the Union2 data set by 
taking into account a global minimization of the five pa- 
rameters. The latter procedure leads to the best fit val- 
ues (h,a,P,Q c ,n x ) = (0.83, 1.08,0.005,0.21,0.75) along 
with Xd o f = 0-98 < 1, indicating that the estimation 
made with the Union2 data is trustworthy and consis- 
tent with the one obtained from the Hubble data set. 

The kinematic analysis based on the behavior of de- 
celeration parameter indicates that at z t = 0.87 the uni- 
verse begins to accelerate, being qo = —0.45 its current 
value within la C.L.. Concerning the equations of state, 
we have found that —1 < oj x (z),oj(z) < 0, ui c = 0, 
and ojf, = a — 1 > along with theirs current values 
ujq = —0.63, u! c q = 0, ui x o = —0.82, and cj&o = 0.05 (cf. 
Fig. 2]) . In addition, we have obtained that the age of the 
universe is 13.17 Gyr very close to the one reporte d by 
WMAP-7year project or WMAP +BAO+H data [Hf. 
So that the MHRDE is free from the cosmic-age problem 
at low redshift (0 < z < 2) in contrast to the Ricci-likc 
HDE, giving rise age- redshift curves below two old stellar 
sources such as the 4 Gyr old galaxy LBDS 53W069 at 
redshift z = 1.43 H3 and the 3.5 Gyr old galaxy LBDS 
53W091 at redshift z = 1.55 (cf. Fig. |SJ. 

We have studied the issue of dark energy at early times 
by taking into account the stringent bounds reported at 
recombination era and/or at BBN and shown that the 
inclusion of a non interacting component makes possi- 
ble that the amount of dark energy at Zt ~ 0(1) be- 
gins to decrease sharply, giving fl x (z = 10 3 ) = 0.038 
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and fl x (z ~ 1100) = 0.032 at recombination era. These 
bounds indicate a good agreement with the forecast of 
Planck and CMBPol experiments [33| as well as with the 
upper bound provided by the constraints on the varia- 
tion in the fine structure constant (35j . As we have al- 
ready mentioned, the presence of dark energy at BBN era 
should not disturb the observed Helium abundance in the 
universe which is regarded as one of the major evidence 
in supporting the big bang theory. We have obtained the 
dark energy density parameters Q x (z = 10 9 ) = 0.0068 
and Vl x (z = 10 12 ) = 0.0021 so the fraction of EDE fulfills 
the bound Q x (z ~ 10 10 ) < 0.21 at BBN [H| (see Fig®. 
Finally, we would like to stress that the values of the cos- 
mological parameters obtained here are consistent with 
four important data set: i— the present-day scenario ob- 
tained with the Hubble data and SNe la data set, ii— the 



recombination bounds for EDE, hi— the BBN data, and 
iv— LSS. 
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